Genetic determinants of plasma protein levels in the Estonian population.

The proteome holds great potential as an intermediate layer between the genome and phenome. Previous protein quantitative trait locus studies have focused mainly on describing the effects of common genetic variations on the proteome. Here, we assessed the impact of the common and rare genetic variations as well as the copy number variants (CNVs) on 326 plasma proteins measured in up to 500 individuals. We identified 184 cis and 94 trans signals for 157 protein traits, which were further fine-mapped to credible sets for 101 cis and 87 trans signals for 151 proteins. Rare genetic variation contributed to the levels of 7 proteins, with 5 cis and 14 trans associations. CNVs were associated with the levels of 11 proteins (7 cis and 5 trans), examples including a 3q12.1 deletion acting as a hub for multiple trans associations; and a CNV overlapping NAIP, a sensor component of the NAIP-NLRC4 inflammasome which is affecting pro-inflammatory cytokine interleukin 18 levels. In summary, this work presents a comprehensive resource of genetic variation affecting the plasma protein levels and provides the interpretation of identified effects.

An original aneuploidy-related gene model for predicting lung adenocarcinoma survival and guiding therapy.

Aneuploidy is a hallmark of cancers, but the role of aneuploidy-related genes in lung adenocarcinoma (LUAD) and their prognostic value remain elusive. Gene expression and copy number variation (CNV) data were enrolled from TCGA and GEO database. Consistency clustering analysis was performed for molecular cluster. Tumor microenvironment was assessed by the xCell and ESTIMATE algorithm. Limma package was used for selecting differentially expressed genes (DEGs). LASSO and stepwise multivariate Cox regression analysis were used to establish an aneuploidy-related riskscore (ARS) signature. GDSC database was conducted to predict drug sensitivity. A nomogram was designed by rms R package. TCGA-LUAD patients were stratified into 3 clusters based on CNV data. The C1 cluster displayed the optimal survival advantage and highest inflammatory infiltration. Based on integrated intersecting DEGs, we constructed a 6-gene ARS model, which showed effective prediction for patient's survival. Drug sensitivity test predicted possible sensitive drugs in two risk groups. Additionally, the nomogram exhibited great predictive clinical treatment benefits. We established a 6-gene aneuploidy-related signature that could effectively predict the survival and therapy for LUAD patients. Additionally, the ARS model and nomogram could offer guidance for the preoperative estimation and postoperative therapy of LUAD.

Expanding deep phenotypic spectrum associated with atypical pathogenic structural variations overlapping 15q11-q13 imprinting region.

BACKGROUND: The 15q11-q13 region is a genetic locus with genes subject to genomic imprinting, significantly influencing neurodevelopment. Genomic imprinting is an epigenetic phenomenon that causes differential gene expression based on the parent of origin. In most diploid organisms, gene expression typically involves an equal contribution from both maternal and paternal alleles, shaping the phenotype. Nevertheless, in mammals, including humans, mice, and marsupials, the functional equivalence of parental alleles is not universally maintained. Notably, during male and female gametogenesis, parental alleles may undergo differential marking or imprinting, thereby modifying gene expression without altering the underlying DNA sequence. Neurodevelopmental disorders, such as Prader-Willi syndrome (PWS) (resulting from the absence of paternally expressed genes in this region), Angelman syndrome (AS) (associated with the absence of the maternally expressed UBE3A gene), and 15q11-q13 duplication syndrome (resulting from the two common forms of duplications-either an extra isodicentric 15 chromosome or an interstitial 15 duplication), are the outcomes of genetic variations in this imprinting region. METHODS: Conducted a genomic study to identify the frequency of pathogenic variants impacting the 15q11-q13 region in an ethnically homogenous population from Bangladesh. Screened all known disorders from the DECIPHER database and identified variant enrichment within this cohort. Using the Horizon analysis platform, performed enrichment analysis, requiring at least >60% overlap between a copy number variation and a disorder breakpoint. Deep clinical phenotyping was carried out through multiple examination sessions to evaluate a range of clinical symptoms. RESULTS: This study included eight individuals with clinically suspected PWS/AS, all previously confirmed through chromosomal microarray analysis, which revealed chromosomal breakpoints within the 15q11-q13 region. Among this cohort, six cases (75%) exhibited variable lengths of deletions, whereas two cases (25%) showed duplications. These included one type 2 duplication, one larger atypical duplication, one shorter type 2 deletion, one larger type 1 deletion, and four cases with atypical deletions. Furthermore, thorough clinical assessments led to the diagnosis of four PWS patients, two AS patients, and two individuals with 15q11-q13 duplication syndrome. CONCLUSION: Our deep phenotypic observations identified a spectrum of clinical features that overlap and are unique to PWS, AS, and Dup15q syndromes. Our findings establish genotype-phenotype correlation for patients impacted by variable structural variations within the 15q11-q13 region.

The prevalence of copy number increase at multiallelic copy number variants associated with cave colonization.

Copy number variation is a common contributor to phenotypic diversity, yet its involvement in ecological adaptation is not easily discerned. Instances of parallelly evolving populations of the same species in a similar environment marked by strong selective pressures present opportunities to study the role of copy number variants (CNVs) in adaptation. By identifying CNVs that repeatedly occur in multiple populations of the derived ecotype and are not (or are rarely) present in the populations of the ancestral ecotype, the association of such CNVs with adaptation to the novel environment can be inferred. We used this paradigm to identify CNVs associated with recurrent adaptation of the Mexican tetra (Astyanax mexicanus) to cave environment. Using a read-depth approach, we detected CNVs from previously re-sequenced genomes of 44 individuals belonging to two ancestral surfaces and three derived cave populations. We identified 102 genes and 292 genomic regions that repeatedly diverge in copy number between the two ecotypes and occupy 0.8% of the reference genome. Functional analysis revealed their association with processes previously recognized to be relevant for adaptation, such as vision, immunity, oxygen consumption, metabolism, and neural function and we propose that these variants have been selected for in the cave or surface waters. The majority of the ecotype-divergent CNVs are multiallelic and display copy number increases in cavefish compared to surface fish. Our findings suggest that multiallelic CNVs - including gene duplications - and divergence in copy number provide a fast route to produce novel phenotypes associated with adaptation to subterranean life.

Biallelic variants identified in 36 Pakistani families and trios with autism spectrum disorder.

With its high rate of consanguineous marriages and diverse ethnic population, little is currently understood about the genetic architecture of autism spectrum disorder (ASD) in Pakistan. Pakistan has a highly ethnically diverse population, yet with a high proportion of endogamous marriages, and is therefore anticipated to be enriched for biallelic disease-relate variants. Here, we attempt to determine the underlying genetic abnormalities causing ASD in thirty-six small simplex or multiplex families from Pakistan. Microarray genotyping followed by homozygosity mapping, copy number variation analysis, and whole exome sequencing were used to identify candidate. Given the high levels of consanguineous marriages among these families, autosomal recessively inherited variants were prioritized, however de novo/dominant and X-linked variants were also identified. The selected variants were validated using Sanger sequencing. Here we report the identification of sixteen rare or novel coding variants in fifteen genes (ARAP1, CDKL5, CSMD2, EFCAB12, EIF3H, GML, NEDD4, PDZD4, POLR3G, SLC35A2, TMEM214, TMEM232, TRANK1, TTC19, and ZNF292) in affected members in eight of the families, including ten homozygous variants in four families (nine missense, one loss of function). Three heterozygous de novo mutations were also identified (in ARAP1, CSMD2, and NEDD4), and variants in known X-linked neurodevelopmental disorder genes CDKL5 and SLC35A2. The current study offers information on the genetic variability associated with ASD in Pakistan, and demonstrates a marked enrichment for biallelic variants over that reported in outbreeding populations. This information will be useful for improving approaches for studying ASD in populations where endogamy is commonly practiced.

Evolving copy number gains promote tumor expansion and bolster mutational diversification.

The timing and fitness effect of somatic copy number alterations (SCNA) in cancer evolution remains poorly understood. Here we present a framework to determine the timing of a clonal SCNA that encompasses multiple gains. This involves calculating the proportion of time from its last gain to the onset of population expansion (lead time) as well as the proportion of time prior to its first gain (initiation time). Our method capitalizes on the observation that a genomic segment, while in a specific copy number (CN) state, accumulates point mutations proportionally to its CN. Analyzing 184 whole genome sequenced samples from 75 patients across five tumor types, we commonly observe late gains following early initiating events, occurring just before the clonal expansion relevant to the sampling. These include gains acquired after genome doubling in more than 60% of cases. Notably, mathematical modeling suggests that late clonal gains may contain final-expansion drivers. Lastly, SCNAs bolster mutational diversification between subpopulations, exacerbating the circle of proliferation and increasing heterogeneity.

Allele-specific transcriptional effects of subclonal copy number alterations enable genotype-phenotype mapping in cancer cells.

Subclonal copy number alterations are a prevalent feature in tumors with high chromosomal instability and result in heterogeneous cancer cell populations with distinct phenotypes. However, the extent to which subclonal copy number alterations contribute to clone-specific phenotypes remains poorly understood. We develop TreeAlign, which computationally integrates independently sampled single-cell DNA and RNA sequencing data from the same cell population. TreeAlign accurately encodes dosage effects from subclonal copy number alterations, the impact of allelic imbalance on allele-specific transcription, and obviates the need to define genotypic clones from a phylogeny a priori, leading to highly granular definitions of clones with distinct expression programs. These improvements enable clone-clone gene expression comparisons with higher resolution and identification of expression programs that are genomically independent. Our approach sets the stage for dissecting the relative contribution of fixed genomic alterations and dynamic epigenetic processes on gene expression programs in cancer.

Dynamics of Mitochondrial DNA Copy Number and Membrane Potential in Mouse Pre-Implantation Embryos: Responses to Diverse Types of Oxidative Stress.

Mitochondria undergo a myriad of changes during pre-implantation embryo development, including shifts in activity levels and mitochondrial DNA (mtDNA) replication. However, how these distinct aspects of mitochondrial function are linked and their responsiveness to diverse stressors is not well understood. Here, we show that mtDNA content increased between 8-cell embryos and the blastocyst stage, with similar copy numbers per cell in the inner cell mass (ICM) and trophectoderm (TE). In contrast, mitochondrial membrane potential (MMP) was higher in TE than ICM. Culture in ambient oxygen (20% O2) altered both aspects of mitochondrial function: the mtDNA copy number was upregulated in ICM, while MMP was diminished in TE. Embryos cultured in 20% O2 also exhibited delayed development kinetics, impaired implantation, and reduced mtDNA levels in E18 fetal liver. A model of oocyte mitochondrial stress using rotenone showed only a modest effect on on-time development and did not alter the mtDNA copy number in ICM; however, following embryo transfer, mtDNA was higher in the fetal heart. Lastly, endogenous mitochondrial dysfunction, induced by maternal age and obesity, altered the blastocyst mtDNA copy number, but not within the ICM. These results demonstrate that mitochondrial activity and mtDNA content exhibit cell-specific changes and are differentially responsive to diverse types of oxidative stress during pre-implantation embryogenesis.

Evaluation of the Copy Number Variants and Single-Nucleotide Polymorphisms of ABCA3 in Newborns with Respiratory Distress Syndrome-A Pilot Study.

Background and Objectives: Respiratory distress syndrome (RDS) in preterm infants commonly occurs due to the immaturity-related deficiency of pulmonary surfactant. Beyond prematurity, various environmental and genetic factors can influence the onset and progression of RDS. This study aimed to analyze three single-nucleotide polymorphisms (SNPs) of the ABCA3 gene to assess the ABCA3 gene as a candidate gene for susceptibility to RDS and overall survival in newborns and to evaluate the utility of MLPA in RDS neonatal patients. Materials and Methods: Three SNPs were chosen and genotyped in a cohort of 304 newborns. Data analysis and statistical tests were employed to examine allele frequencies, haplotypes, and measures of pairwise linkage disequilibrium. Results: There was no observed haplotype association with SNPs rs13332514 (c.1059G>A) and rs170447 (c.1741+33T>C) among newborns, both with and without RDS (p > 0.05). The minor C allele frequency of the ABCA3 rs323043 (c.1755G>C) SNP showed a significant increase in preterm infants with RDS. MLPA results indicated that the predominant findings were normal, revealing no CNVs in the genes ABCA3 and SFTPC that were investigated in our patients. Conclusions: The presence of the variant C allele in the rs323043 (c.1755G>C) SNP may be a risk factor for RDS in premature newborns.

Use of Enzymatically Converted Cell-Free DNA (cfDNA) Data for Copy Number Variation-Linked Fragmentation Analysis Allows for Early Colorectal Cancer Detection.

The use of non-invasive liquid biopsy-based cell-free DNA (cfDNA) analysis is an emerging method of cancer detection and intervention. Different analytical methodologies are used to investigate cfDNA characteristics, resulting in costly and long analysis processes needed for combining different data. This study investigates the possibility of using cfDNA data converted for methylation analysis for combining the cfDNA fragment size with copy number variation (CNV) in the context of early colorectal cancer detection. Specifically, we focused on comparing enzymatically and bisulfite-converted data for evaluating cfDNA fragments belonging to chromosome 18. Chromosome 18 is often reported to be deleted in colorectal cancer. We used counts of short and medium cfDNA fragments of chromosome 18 and trained a linear model (LDA) on a set of 2959 regions to predict early-stage (I-IIA) colorectal cancer on an independent test set. In total, 87.5% sensitivity and 92% specificity were obtained on the enzymatically converted libraries. Repeating the same workflow on bisulfite-converted data yielded lower accuracy results with 58.3% sensitivity, implying that enzymatic conversion preserves the cancer fragmentation footprint in whole genome data better than bisulfite conversion. These results could serve as a promising new avenue for the early detection of colorectal cancer using fragmentation and methylation approaches on the same datasets.

Copy number variations within fibroblast growth factor 13 gene influence growth traits and alternative splicing in cattle.

Previous researches revealed a copy number variation (CNV) region in the bovine fibroblast growth factor 13 (FGF13) gene. However, its effects remain unknown. This study detected the various copy number types in seven Chinese cattle breeds and analysed their population genetic characteristics and effects on growth traits and transcription levels. Copy number Loss was more frequent in Caoyuan Red cattle and Xianan cattle than in the other breeds. Association analysis between CNV and growth traits of Qinchuan indicated that the CNV was significantly related to chest depth, hip width and hucklebone width (P < 0.05). Additionally, the growth traits of individuals with copy number Loss were significantly inferior to those with copy number Gain or Median (P < 0.05). Besides, we found two splicing isoforms, AS1 and AS2, in FGF13 gene, which resulted from alternative 5' splicing sites of intron 1. These isoforms showed varied expression levels in various tissues. Moreover, CNV was significantly and negatively associated with the mRNA expression of AS1 (r = -0.525, P < 0.05). The CNVs in bovine FGF13 gene negatively regulated growth traits and gene transcription. These observations provide new insights into bovine FGF13 gene, delivering potentially useful information for future Chinese cattle breeding programs.

Copy number variation sequencing for the products of conception: What is the optimal testing strategy.

BACKGROUND: Copy number variation sequencing (CNV-seq) is crucial in prenatal diagnosis, but its limitations in detecting polyploidy, maternal cell contamination (MCC), and uniparental disomy (UPD) restrict its application in the analysis of products of conception (POCs). This study aimed to investigate an optimal genetic testing strategy for POCs in the era of CNV-seq. METHODS: CNV-seq and quantitative fluorescent polymerase chain reaction (QF-PCR) were performed in all 4,211 spontaneous miscarriage cases. Different testing strategies were compared and the optimal testing strategies were proposed. RESULTS: Of the 4,211 cases, 2561 (60.82%) exhibited clinically significant chromosomal abnormalities. CNV-seq alone, without QF-PCR, might misdiagnose 311 (7.39%) cases, including 278 polyploidy, 13 UPD, and 20 MCC. In 20 MCC cases identified by QF-PCR, CNV-seq successfully pinpointed the cause of miscarriage in 13 cases. Furthermore, in cases where QF-PCR suggested polyploidy, CNV-seq improved the diagnostic accuracy in 54 (1.28%) hypo/hypertriploidy cases. After comparing four different strategies, the sequential approach (initiating with CNV-seq followed by QF-PCR if necessary) emerged as advantageous, reducing approximately 70% of the cost associated with QF-PCR while maintaining result accuracy. CONCLUSIONS: We propose an initial CNV-seq followed by QF-PCR if needed-an efficient and cost-effective strategy for the genetic analysis of POCs.

Psychosis spectrum symptoms among individuals with schizophrenia-associated copy number variants and evidence of cerebellar correlates of symptom severity.

The 3q29 deletion (3q29Del) is a copy number variant (CNV) with one of the highest effect sizes for psychosis-risk (>40-fold). Systematic research offers avenues for elucidating mechanism; however, compared to CNVs like 22q11.2Del, 3q29Del remains understudied. Emerging findings indicate that posterior fossa abnormalities are common among carriers, but their clinical relevance is unclear. We report the first in-depth evaluation of psychotic symptoms in participants with 3q29Del (N=23), using the Structured Interview for Psychosis-Risk Syndromes, and compare this profile to 22q11.2Del (N=31) and healthy controls (N=279). We also explore correlations between psychotic symptoms and posterior fossa abnormalities. Cumulatively, 48% of the 3q29Del sample exhibited a psychotic disorder or attenuated positive symptoms, with a subset meeting criteria for clinical high-risk. 3q29Del had more severe ratings than controls on all domains and only exhibited less severe ratings than 22q11.2Del in negative symptoms; ratings demonstrated select sex differences but no domain-wise correlations with IQ. An inverse relationship was identified between positive symptoms and cerebellar cortex volume in 3q29Del, documenting the first clinically-relevant neuroanatomical connection in this syndrome. Our findings characterize the profile of psychotic symptoms in the largest 3q29Del sample reported to date, contrast with another high-impact CNV, and highlight cerebellar involvement in psychosis-risk.

Digital droplet PCR versus quantitative PCR for lipoprotein (a) kringle IV type 2 repeat polymorphism genetic characterization.

BACKGROUND: Lipoprotein(a) [Lp(a)] level variability, related to atherothrombotic risk increase, is mainly attributed to LPA gene, encoding apolipoprotein(a), with kringle IV type 2 (KIV2) copy number variation (CNV) acting as the primary genetic determinant. Genetic characterization of Lp(a) is in continuous growth; nevertheless, the peculiar structural characteristics of this variant constitute a significant challenge to the development of effective detection methods. The aim of the study was to compare quantitative real-time PCR (qPCR) and digital droplet PCR (ddPCR) in the evaluation of KIV2 repeat polymorphism. METHODS: We analysed 100 subjects tested for cardiovascular risk in which Lp(a) plasma levels were assessed. RESULTS: Correlation analysis between CNV values obtained with the two methods was slightly significant (R = 0.413, p = 0.00002), because of the wider data dispersion in qPCR compared with ddPCR. Internal controls C1, C2 and C3 measurements throughout different experimental sessions revealed the superior stability of ddPCR, which was supported by a reduced intra/inter-assay coefficient of variation determined in this method compared to qPCR. A significant inverse correlation between Lp(a) levels and CNV values was confirmed for both techniques, but it was higher when evaluated by ddPCR than qPCR (R = -0.393, p = 0.000053 vs R = -0.220, p = 0.028, respectively). When dividing subjects into two groups according to 500 mg/L Lp(a) cut-off value, a significantly lower number of KIV2 repeats emerged among subjects with greater Lp(a) levels, with stronger evidence in ddPCR than in qPCR (p = 0.000013 and p = 0.001, respectively). CONCLUSIONS: Data obtained support a better performance of ddPCR in the evaluation of KIV2 repeat polymorphism.

Detection of elusive DNA copy-number variations in hereditary disease and cancer through the use of noncoding and off-target sequencing reads.

Copy-number variants (CNVs) play a substantial role in the molecular pathogenesis of hereditary disease and cancer, as well as in normal human interindividual variation. However, they are still rather difficult to identify in mainstream sequencing projects, especially involving exome sequencing, because they often occur in DNA regions that are not targeted for analysis. To overcome this problem, we developed OFF-PEAK, a user-friendly CNV detection tool that builds on a denoising approach and the use of "off-target" DNA reads, which are usually discarded by sequencing pipelines. We benchmarked OFF-PEAK on data from targeted sequencing of 96 cancer samples, as well as 130 exomes of individuals with inherited retinal disease from three different populations. For both sets of data, OFF-PEAK demonstrated excellent performance (>95% sensitivity and >80% specificity vs. experimental validation) in detecting CNVs from in silico data alone, indicating its immediate applicability to molecular diagnosis and genetic research.

Protein-altering variants at copy number-variable regions influence diverse human phenotypes.

Copy number variants (CNVs) are among the largest genetic variants, yet CNVs have not been effectively ascertained in most genetic association studies. Here we ascertained protein-altering CNVs from UK Biobank whole-exome sequencing data (n = 468,570) using haplotype-informed methods capable of detecting subexonic CNVs and variation within segmental duplications. Incorporating CNVs into analyses of rare variants predicted to cause gene loss of function (LOF) identified 100 associations of predicted LOF variants with 41 quantitative traits. A low-frequency partial deletion of RGL3 exon 6 conferred one of the strongest protective effects of gene LOF on hypertension risk (odds ratio = 0.86 (0.82-0.90)). Protein-coding variation in rapidly evolving gene families within segmental duplications-previously invisible to most analysis methods-generated some of the human genome's largest contributions to variation in type 2 diabetes risk, chronotype and blood cell traits. These results illustrate the potential for new genetic insights from genomic variation that has escaped large-scale analysis to date.

Comprehensive whole-genome sequence analyses provide insights into the genomic architecture of cerebral palsy.

We performed whole-genome sequencing (WGS) in 327 children with cerebral palsy (CP) and their biological parents. We classified 37 of 327 (11.3%) children as having pathogenic/likely pathogenic (P/LP) variants and 58 of 327 (17.7%) as having variants of uncertain significance. Multiple classes of P/LP variants included single-nucleotide variants (SNVs)/indels (6.7%), copy number variations (3.4%) and mitochondrial mutations (1.5%). The COL4A1 gene had the most P/LP SNVs. We also analyzed two pediatric control cohorts (n = 203 trios and n = 89 sib-pair families) to provide a baseline for de novo mutation rates and genetic burden analyses, the latter of which demonstrated associations between de novo deleterious variants and genes related to the nervous system. An enrichment analysis revealed previously undescribed plausible candidate CP genes (SMOC1, KDM5B, BCL11A and CYP51A1). A multifactorial CP risk profile and substantial presence of P/LP variants combine to support WGS in the diagnostic work-up across all CP and related phenotypes.

Effects of gene dosage and development on subcortical nuclei volumes in individuals with 22q11.2 copy number variations.

The 22q11.2 locus contains genes critical for brain development. Reciprocal Copy Number Variations (CNVs) at this locus impact risk for neurodevelopmental and psychiatric disorders. Both 22q11.2 deletions (22qDel) and duplications (22qDup) are associated with autism, but 22qDel uniquely elevates schizophrenia risk. Understanding brain phenotypes associated with these highly penetrant CNVs can provide insights into genetic pathways underlying neuropsychiatric disorders. Human neuroimaging and animal models indicate subcortical brain alterations in 22qDel, yet little is known about developmental differences across specific nuclei between reciprocal 22q11.2 CNV carriers and typically developing (TD) controls. We conducted a longitudinal MRI study in a total of 385 scans from 22qDel (n = 96, scans = 191, 53.1% female), 22qDup (n = 37, scans = 64, 45.9% female), and TD controls (n = 80, scans = 130, 51.2% female), across a wide age range (5.5-49.5 years). Volumes of the thalamus, hippocampus, amygdala, and anatomical subregions were estimated using FreeSurfer, and the linear effects of 22q11.2 gene dosage and non-linear effects of age were characterized with generalized additive mixed models (GAMMs). Positive gene dosage effects (volume increasing with copy number) were observed for total intracranial and whole hippocampus volumes, but not whole thalamus or amygdala volumes. Several amygdala subregions exhibited similar positive effects, with bi-directional effects found across thalamic nuclei. Distinct age-related trajectories were observed across the three groups. Notably, both 22qDel and 22qDup carriers exhibited flattened development of hippocampal CA2/3 subfields relative to TD controls. This study provides novel insights into the impact of 22q11.2 CNVs on subcortical brain structures and their developmental trajectories.

Somatic CDKN2A copy number variations are associated with the prognosis of esophageal squamous cell dysplasia.

BACKGROUND: Somatic copy number variations (SCNVs) in the CDKN2A gene are among the most frequent events in the dysplasia-carcinoma sequence of esophageal squamous cell carcinoma. However, whether CDKN2A SCNVs are useful biomarkers for the risk stratification and management of patients with esophageal squamous cell dysplasia (ESCdys) is unknown. This study aimed to investigate the characteristics and prognostic value of CDKN2A SCNVs in patients with mild or moderate (m/M) ESCdys. METHODS: This study conducted a prospective multicenter study of 205 patients with a baseline diagnosis of m/M ESCdys in five high-risk regions of China (Ci County, Hebei Province; Yanting, Sichuan Province; Linzhou, Henan Province; Yangzhong, Jiangsu Province; and Feicheng, Shandong Province) from 2005 to 2019. Genomic DNA was extracted from paraffin biopsy samples and paired peripheral white blood cells from patients, and a quantitative polymerase chain reaction assay, P16-Light, was used to detect CDKN2A copy number. The cumulative regression and progression rates of ESCdys were evaluated using competing risk models. RESULTS: A total of 205 patients with baseline m/M ESCdys were enrolled. The proportion of ESCdys regression was significantly lower in the CDKN2A deletion cohort than in the diploid and amplification cohorts (18.8% [13/69] vs. 35.0% [28/80] vs. 51.8% [29/56], P  <0.001). In the univariable competing risk analysis, the cumulative regression rate was statistically significantly lower ( P = 0.008), while the cumulative progression rate was higher ( P = 0.017) in ESCdys patients with CDKN2A deletion than in those without CDKN2A deletion. CDKN2A deletion was also an independent predictor of prognosis in ESCdys ( P = 0.004) in the multivariable analysis. CONCLUSION: The results indicated that CDKN2A SCNVs are associated with the prognosis of ESCdys and may serve as potential biomarkers for risk stratification.

A method to correct for local alterations in DNA copy number that bias functional genomics assays applied to antibiotic-treated bacteria.

Functional genomics techniques, such as transposon insertion sequencing and RNA-sequencing, are key to studying relative differences in bacterial mutant fitness or gene expression under selective conditions. However, certain stress conditions, mutations, or antibiotics can directly interfere with DNA synthesis, resulting in systematic changes in local DNA copy numbers along the chromosome. This can lead to artifacts in sequencing-based functional genomics data when comparing antibiotic treatment to an unstressed control. Further, relative differences in gene-wise read counts may result from alterations in chromosomal replication dynamics, rather than selection or direct gene regulation. We term this artifact "chromosomal location bias" and implement a principled statistical approach to correct it by calculating local normalization factors along the chromosome. These normalization factors are then directly incorporated into statistical analyses using standard RNA-sequencing analysis methods without modifying the read counts themselves, preserving important information about the mean-variance relationship in the data. We illustrate the utility of this approach by generating and analyzing a ciprofloxacin-treated transposon insertion sequencing data set in Escherichia coli as a case study. We show that ciprofloxacin treatment generates chromosomal location bias in the resulting data, and we further demonstrate that failing to correct for this bias leads to false predictions of mutant drug sensitivity as measured by minimum inhibitory concentrations. We have developed an R package and user-friendly graphical Shiny application, ChromoCorrect, that detects and corrects for chromosomal bias in read count data, enabling the application of functional genomics technologies to the study of antibiotic stress.IMPORTANCEAltered gene dosage due to changes in DNA replication has been observed under a variety of stresses with a variety of experimental techniques. However, the implications of changes in gene dosage for sequencing-based functional genomics assays are rarely considered. We present a statistically principled approach to correcting for the effect of changes in gene dosage, enabling testing for differences in the fitness effects or regulation of individual genes in the presence of confounding differences in DNA copy number. We show that failing to correct for these effects can lead to incorrect predictions of resistance phenotype when applying functional genomics assays to investigate antibiotic stress, and we provide a user-friendly application to detect and correct for changes in DNA copy number.